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ABSTRACT 

Certain  elastic  solids  when  subjected  to  sufficiently  high  loads  abruptly 
change  their  mechanical  properties  and  yet  continue  to  respond  elastically  to 
further  loading.  In  one  dimension  such  mechanically  induced  elastic  phase 
transitions  may  be  due  to  a  non-monotonic  stress-strain  relationship.  This 
appears  to  be  particularly  true  for  certain  mineral  crystals,  such  as  calcite. 

This  work  considers  a  one-dimensional  dynamical  problem  for  a  special 
material.  The  problem  reduces  to  determining  the  location  of  the  internal 
moving  boundary  separating  distinct  elastic  phases.  This  phase  boundary  is 
similar  to  a  gas  dynamical  shock  wave.  For  the  problem  considered  here,  this 
phase  boundary  is  shown  to  be  governed  by  a  functional  equation  of  the  form 

*<4>(t))  +  F{$(t) )  +  t  -  0 

for  the  unknown  $(t),  where  F( }  is  a  known  function  involving  the  boundary 
conditions.  Hie  unusual  equation  is  derived  by  considering  the  effect  of 
acoustic  waves  repeatedly  reflecting  between  the  phase  boundary  and  the 
external  boundary.  Hie  equation  is  shown  to  possess  a  unique  solution  and  is 
treated  asymptotically  to  determine  the  large-time  behavior  of  the  phase 
boundary. 
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SIGNIFICANCE  AND  EXPLANATION 


Certain  elastic  solids  when  subjected  to  sufficiently  high  loads  abruptly 
change  their  mechanical  properties  and  yet  continue  to  respond  elastically  to 
further  loading.  In  one  dimension  such  mechanically  induced  el i., tic  phase 
transitions  may  be  due  to  a  non-monotonic  stress-strain  curve.  This  -work 
investigates  the  cumulative  reflection  of  acoustic  waves  between  the  external 
boundary  of  the  solid  and  the  internal  moving  boundary  separating  distinct 
elastic  phases.  This  latter  phase  boundary  is  similar  to  a  gas  dynamical 
shock  wave.  For  the  material  introduced  in  this  work,  a  functional  equation 
governing  the  trajectory  of  a  phase  boundary  is  derived  and  shown  to  have  a 
unique  solution.  This  equation  is  treated  asymptotically  to  determine  the 
large  time  behavior  of  the  phase  boundary. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC  and  not  with  the  author  of  this  report. 


A  FUNCTIONAL  EQUATION  GOVERNING  MOVING  PHASE  BOUNDARIES  IN  AN  ELASTIC  BAR 

Thomas  J.  Pence 

1 .  Introduction 

Certain  elastic  solids  when  subjected  to  sufficiently  high  loads  will  abruptly  change 
their  mechanical  properties  and  yet  continue  to  respond  elastically  to  further  loading. 

An  example  of  this  type  of  behavior  is  found  in  the  mineral  calcite,  which  suffers  a 
change  in  its  crystal  structure  whenever  it  is  sufficiently  compressed  [1].  James  [2], 

[3]  cites  other  examples  in  materials  ranging  from  metals  to  natural  rubbers  and  various 
polymers.  The  possibility  that  a  material  may  cease  to  behave  in  an  expected  manner  has 
an  impact  on  a  host  of  geophysical  and  engineering  problems.  It  is  also  interesting  to 
note  that  metals  which  undergo  a  phase  transition  due  to  a  change  in  temperature  are  being 
promoted  as  a  means  for  performing  specific  engineering  tasks  [4] . 

The  study  of  such  phenomena  within  the  theory  of  elasticity  centers  around  material 
models  involving  non-convex  strain  energies  (in  addition  to  [2],  [3]  see  Ericksen  [5], 
Knowles  &  Sternberg  [6],  [7],  [8],  Abeyaratne  [9],  [10]  and  Pence  [11]).  In  one 
dimension,  these  models  are  equivalent  to  theories  of  materials  with  a  non-monotonic 
stress-strain  relation.  This  type  of  model  is  not  new;  the  theory  of  a  Van  der  Waals 
fluid  is  based  upon  a  non-monotonic  constitutive  relation.  Indeed,  phase  boundaries  in 
certain  solids  may  be  kinematically  similar  to  the  liquid-vapor  phase  boundary  in  a  Van 
der  Waals  fluid.  Admissibility  criteria  for  the  later  type  of  phase  boundary  have  been 
studied  by  Serrin  [12],  Slemrod  [13],  and  Slemrod  &  Hagan  [14]. 

This  paper  considers  the  problem  of  determining  the  motion  of  the  front  associated 
with  a  phase  transition  for  a  model  elastic  material  in  one  dimension,  we  imagine  a  bar 
in  which  transverse  displacements  are  absent  and  we  let  u(x,t)  be  the  longitudinal 

displacement  in  a  Lagrangian  frame.  We  also  let  s,  e  ■  -|^  and  v  *  -|^  denote 
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respectively  the  stress,  strain  and  velocity  in  the  bar.  The  material  is  assumed  to  obey 
a  stress-strain  relation  like  that  depicted  in  figure  1. 


a(e) 


Figure  1. 


The  function  o(e)  is  taken  to  be  of  unlimited  smoothness,  concave  dam  on  the  first 
ascending  branch,  and  is  assumed  to  approach  “  as  e  •*■  “>.  These  restrictions  are  not 
essential  to  the  general  theory;  what  is  essential  is  the  existence  of  two  distinct 
ascending  branches.  It  is  these  two  branches  which  are  associated  with  distinct  material 
phases.  We  consider  a  tensile  load  g(t)  applied  at  x  ■  0  to  a  bar  which  is  initially 
undeformed  and  at  rest,  and  which  is  also  taken  to  be  semi-infinite  in  order  to  avoid 
having  to  treat  waves  reflected  from  a  second  fixed  boundary.  The  equation  of  motion  is 


(1.1) 

It  is  to  be  solved  subject  to 


32u  3a 


(1.2)  o(e(0,t) )  -  g(t),  (t  >  0), 

(1.3)  c(x,0)  -  0,  v(x, 0)  =0,  (x  >  0)  . 
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As  discussed  in  [11],  solutions  will  include  discontinuities  in  a,  e,  v  in  the  event 
that  g(t)  >  Og.  If  x  =  s(t)  is  the  trajectory  in  the  (x,t)-plane  of  such  a 
discontinuity,  the  jump  in  the  dynamical  fields  across  x  =  s(t)  are  to  be  restricted  by 
the  familiar  conditions 

(1.4,  If  i  «  I  ♦  I  V  I  -  0,  If  1  V  1  +  |  o  |  -  0  . 

Phase  boundaries  are  those  particular  discontinuity  fronts  that  separate  strains 
associated  with  the  two  different  ascending  branches  of  fig.  1. 

The  model  is  explained  more  fully  in  section  2,  where  we  also  consider  the  special 
problem  of  an  impulsive  load  (g(t)  a  step  function).  This  problem  is  studied  in  order  to 
touch  upon  some  questions  concerning  uniqueness  as  well  as  in  order  to  discuss  tun  unusual 
-  and  physically  desirable  -  feature  of  certain  mathematical  solutions:  phase  boundaries 
which  travel  at  lower  speeds  than  the  sound  speeds  of  the  material  in  both  adjacent 
phases.  In  section  3,  we  return  to  a  more  general  loading  program,  and  introduce  the 
notion  of  an  energy  confining  phase  boundary.  In  section  4  we  restrict  attention  to 
materials  for  which  the  second  ascending  branch  is  linear.  This  gives  rise  to  a  problem 
involving  a  linear  wave  equation  in  a  region  of  the  (x.t)-plane  which  is  partially  bounded 
by  the  (as  yet  unknown)  location  of  a  phase  boundary.  This  problem  is  reducible  in 
certain  cases  to  a  functional  equation  of  the  form 
(1.5)  <t>($U))  "  2$(z)  +  z  «  Q(*(z>)  , 

where  $,  the  function  to  be  determined,  is  subject  to  some  additional  conditions  and 
Q  is  a  known  monotonically  increasing  function  which  incorporates  the  function  g(t). 

The  unusual  equation  (1.5)  is  derived  by  considering  the  cumulative  effect  of  acoustic 
waves  repeatedly  reflecting  between  the  moving  phase  boundary  and  the  fixed  external 
boundary.  In  section  5,  we  prove  the  existence  of  a  unique  solution  $  by  means  of  an 
iterative  scheme.  A  more  delicate  argument  subsequently  guarantees  the  differentiability 
of  this  solution.  In  the  final  section,  eqn.  (1.5)  is  treated  asymptotically  to  determine 
the  large-time  behavior  of  the  phase  boundary.  The  resulting  asymptotic  solution  agrees 
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with  the  particular  solution  of  the  impulsive  load  problem  presented  in  section  2  in  the 
limit  t  ♦  ",  x/t  fixed. 

Equations  in  which  a  function  to  be  determined  appears  in  its  own  argument  arise  in 
other  branches  of  mathematics.  Specifically,  (1.5)  is  a  special  case  of  the  equation  of 
invariant  curves.  Thus  the  theoretical  results  presented  in  section  5  extend  some  pre¬ 
vious  results  in  that  theory  (see  Chapter  XIV  of  [15]).  In  the  field  of  continuum  me¬ 
chanics,  problems  involving  wave  propagation  can  lead  to  equations  involving  a  delay  in 
their  argument.  Although  this  work  is  devoted  to  solid  mechanics,  a  set  of  functional 
equations  have  been  considered  by  Seymour  &  Mortel  [16]  in  connection  with  the  study  of 
oscillations  occurring  in  an  inviscid  gas  confined  in  a  closed  tube.  In  fact,  their  eqn. 
1.  can  be  manipulated  into  the  form  of  (1.5),  although  the  properties  of  Q  (H  in  [16]) 
differ.  In  both  this  work  and  [16],  the  appearance  of  a  function  in  its  own  argument  is 
due  to  the  unknown  time  delay  for  the  return  of  a  reflected  signal.  It  is  intriguing  that 
in  [16]  this  delay  is  due  to  an  amplitude  dependent  acoustic  speed,  whereas  in  this  work 
the  sound  speeds  are  known  prior  to  the  derivation  of  (1.5).  In  the  problem  considered 
here,  the  unknown  delay  stems  from  the  unknown  location  of  the  reflecting  phase  boundary. 


2.  Solutions  for  the  Impulsive  Load  Problem 


The  sound  speed  of  the  material  is  given  by  the  value  /o‘(c>,  where  the  '  symbol 

is  the  usual  notation  for  derivative •  Shock  waves  arise  naturally  in  the  solution  of 

(1.1-4)  due  to  the  intersection  of  characteristic  curves  in  the  (x,t)-plane.  In  contrast. 

phase  boundaries  occur  in  this  theory  because  no  single  branch  of  the  c(e) -curve  can 

accommodate  stress  values  both  less  than  a  and  greater  than  o..  As  discussed  in  [3], 

a  p 

phase  boundaries  differ  from  the  aforementioned  shock  waves  in  that  they  may  travel  at  a 
much  lower  velocity  than  the  sound  speeds  of  the  material  in  each  adjacent  phase. 

In  an  equilibrium  setting,  configurations  involving  strains  on  the  descending  branch 
of  the  o(e) -curve  are  extremely  unstable  [2].  For  the  dynamical  problem,  the  sound 
speeds  are  no  longer  real  on  this  branch.  This  prompts  us  to  seek  solutions  to  (1.1-4) 
which  avoid  this  branch  altogether.  Thus  every  time-interval  during  which  g(t) 
traverses  the  range  must  generate  at  least  one  jump  in  e(0,t)  between  the  two 

ascending  branches.  Let  us  suppose  that  g(t)  increases  monotonically  from  0  to  a 
value  greater  than  Og  and  discuss  solutions  with  only  one  such  jump,  say  first  occurring 
at  t  -  b**  Then  a  <  g(t*)  S  o.  <  o.  where  o_  is  the  stress  at  which  the  change  of 

Cl  J  p  u 

phase  first  occurs. 

A  criterion  for  selecting  Oj  remains  elusive.  A  well-known  stability  argument 
yields  the  value  associated  with  the  Maxwell  line  [3] .  Recent  work  in  this  vein  has 
addressed  hew  strsin  rates  and  strain  gradients  can  affect  o,  (see  [12),  113),  [14)  for 

w 

the  corresponding  problem  in  a  Van  der  Waals  gas),  as  well  as  how  a  may  be  affected  by 

«J 

a  phase  boundary  surface  energy  [17].  This  paper  addresses  the  simplier  and  more 
immediate  problem  of  constructing  solutions  to  (1.1-4)  for  a  given  value  of  Thus  we 

append  the  following  restrictions  upon  the  strain  at  x  •  0: 

(2.1)  0  <  e(0,t)  <  tj  or  e(0,t)  >7Jt  (t  >  0)  , 

where  t,  and  e,  denote  respectively  the  first  -  and  second  -  phase  values  of  strain 
J  J 

associated  with  o  via 
J 
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(2.2) 


aUj)  -  o(e j)  -  Oj,  a,  <  Ej  <  B,,  “2  *  eJ  *  62  ' 

One  might  expect  condition  (2.1)  to  render  solutions  of  (1.1-4)  unique,  this,  however,  is 
not  necessarily  true.  Consider  first  the  impulsive-load  problem 


(2.3) 


0, 


0  , 


g(t) 


0.  >  0j,  t  >  0  . 


Following  James  [3],  we  seek  a  similarity  solution  of  the  form 

(2.4)  £(x,t)  »  £(\),  v(x,t)  -  v(A),  X  ■  x/t  , 

where  £(•),  v(*)  are  piecewise  differentiable  functions  defined  on  0  <  X  <  “>.  if 
either  function  is  discontinuous  at  X  «*  X^,  then  x  =  X^t  is  either  a  shock  or  phase 
boundary.  Since  every  ray  x  *  Xt  meets  at  the  origin,  we  shall  in  this  section 
interpret  (2.1)  as  requiring  e(A)  e  (0,£j]  IJ  (£^,“1 .  Equations  for  £  and  v  follow  by 
introducing  (2.4)  into  (1.1).  One  finds  that  £(X)  is  either  constant  or  obeys 

(2.5)  o' (e(A) )  -  X2  . 

The  function  v(A)  satisfies 


(2.6) 


3v  .  3e 
3X  *  'X  IX  ’ 


This  last  equation,  in  conjunction  with  (1.3-4)  indicate  that  v  follows  from  £  as 


£*(  A ) 

X,  <  x  <  -  , 

-  / 

/o' ( 8) 

ds. 

(2.7) 

v(  A )  - 

0 

£  ( A ) 

A  <  A  <  A 

n+1  n 

(n  *  1,  ...,  N  —  1)  , 

v(  X  -) 

'  I 

/C  (s)  ds, 

X 

E(X.  ) 

1 

0  <  x  <  xN  , 

where 

Al  >  X2 

i  ...  >  xN 

are  the 

values  of  X 

at  which  £  is 

discontinuous  and 

(2.8) 

v(XA-)  = 

v(A .+)  +  A 

i  i  1 

[c(Xi+)  -  e ( X±- ) ] 

• 

Thus  it 

merely 

remains  to 

find  an 

e(A)  which: 

(a)  is  either  constant  or  obeys  ( 2.5)  on  some  partitioning  of  0  <  A  <  “  into  open 
intervals  , 
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(b) 


(2.9) 


yields  the  correct  boundary  values  for  e(x,t)  via  the  relations/ 
e(“>)  =  0,  £ ( 0 )  =toc  where  o( e^)  =  ( 


(c)  at  points  of  discontinuity  X^  obeys 

(2.10)  X*  [  e(X.)  J  =  [  o(e(Xi))  ]  , 

and 

(d)  has  range  [0,e._]  'c  [£_,"). 

J  J 

Notice  that  (2.8),  (2.10)  are  necessary  and  sufficient  to  guarantee  the  discontinuity 
conditions  (1.4).  We  now  present  a  method  for  generating  solutions  e(X)  obeying  (a)  - 
(d).  The  constriction  makes  use  of  an  upper  concave  envelope  of  portions  of  the  function 
o(e) . 

A  function  h(x)  defined  on  a  <  x  4  b  is  concave  if  h(Xx  +  (1  -  X)y)  >  Xh(x)  + 

(1  -  X)h(y)  for  all  x,y  e  (a,bj  and  all  X  e  [0,1).  The  upper  concave  envelope  k(x) 
of  a  function  k(x)  is  the  smallest  concave  function  obeying  k(x)  >  k(x).  Thus  at  each 
fixed  x0  either  k(xQ)  =  or  vxo'  k^xo^  Ires  on  a  line  segment  bridging  two 

points  on  the  graph  of  k.  Figure  2  depicts  0(e),  the  upper  concave  envelope  of  0(e) 
in  the  domain  [0,e  ]  for  a  value  e  as  given  in  (2.9). 

CO  OO 


C(c.) 


l- PI  er  ecncavo  envelope-  of  ’  ( f )  given  in  figure  1,  for 
a  rpeoifi:;  vaJut.- 


Figure  2. 


Since  the  derivative  o'(£)  decreases  continuously  from  o'(0)  to  u'  (e^)  the  equation 


(2.11) 


X“  =  a'(e)  , 


has  at  least  one  solution  £  for  each  value  X  e  (0)].  The  value  £  is 

unique  if  £  does  not  lie  on  a  linear  portion  of  o(e).  Let  X^  be  the  positive  square 
root  of  the  slope  of  the  iLtl  linear  portion  of  o(e).  Then  X1  >  X^  >  •••  >  X^  and 
(2.11)  has  an  interval  of  solutions  £  for  each  X  =  X^ .  Thus  (2.11)  implicitly  defines 


a  decreasing  function  £(X)  on 


/o'  (e_) 


<  X  <  /o'(0)  which  is  discontinuous  at  each 


X  *  X^.  By  c  nstruction  e  obeys  (2.5)  between  these  discontinuities.  At  the  same 

2 

time,  (2.10)  holds  for  all  X  =  Xi  on  account  of  X^  being  the  slope  of  a  secant  line  to 

the  curve  o(e).  Since  e(/j*(eai))  =  and  c[/o'  (0) )  =  0,  we  can  generate  a  function 
obeying  (a)  -  (c)  by  extending  e(X)  to  the  complete  domain  X  >  0  via 


(2.12) 


£(  X )  -  0. 


X  >  /  o' (0)  , 


£(X)  *  e^,  0  <  X  <  ✓  o'(£a)  . 

It  remains  to  determine  whether  or  not  (d)  holds.  Before  doing  so,  however,  we  note  that 
the  first  line  segment  of  o(e)  traverses  an  e-interval  including  [ B 1 , B2 1  by  virtue  of 
o"(e)  <0  on  0  <  e  4  0^.  Thus  X1  is  the  phase  boundary  speed.  The  next  N  -  2  shock 
speeds  X  are 


(2.13) 


/o'(e(X  -))  =>  X  =  /o ' ( e( X, +) )  ,  (i  »  2 
i  i  i 


,  ...,  N  -  1)  , 


since  linear  portions  of  0(e)  strictly  interior  to  ((>,£,,)  must  join  smoothly  to  points 
on  0(e).  Thus  the  shocks  x  *  X^t  (i  =  2,  ...,  N  -  1)  are  transonic  with  respect  to 
the  material  flanking  it  on  both  sides.  If  N  >  1,  the  phase  boundary  speed  X1  obeys 


(2.14) 


/o’(e(Xx-))  =  >  /o'(e(X  +))  ,  (N  >  1)  , 


where  the  inequality  is  strict  only  if  e(X^+)  =  0.  in  such  a  case  the  phase  boundary  is 
supersonic  with  respect  to  the  material  in  front.  Similarly  the  fastest  shock  x  =  XNt 
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(2.15) 


>  XN  =  /o'(eUN+)) 


(N  >  1)  . 


Finally,  if  the  phase  boundary  is  the  only  discontinuity  front,  then 
(2.16)  /o’tcU.,-))  >  X1  >  /a 1  (e  ( X ,+  ) }  ,  (N  =  1)  . 


Turning  now  to  the  question  of  whether  or  not  (d)  holds,  we  note  that  nothing  in  the 
above  construction  prevents  e(X)  from  taking  values  in  the  range  t e  , B - ] .  We  can, 
however,  enforce  (d),  while  maintaining  (a)  -  (c),  by  repeating  the  same  construction,  not 
to  the  upper  concave  envelope  of  o(e)  on  0  <  t  <  em,  but  rather  to  the  upper  concave 
envelope  on  0  <  e  <  of  the  "deleted"  function 

^o(e),  0  <  e  <  e  ,  e  <  e  <  ett, 

(2.17)  o(e  ie)  » 

^  not  defined*,  <  e  <  . 

The  program  is  now  to  define  c ( X )  in  a  fashion  analagous  to  (2.11)  by  means  of  the 
equation 

(2.16)  X2  =  o2  ( e e"( X ) )  , 

where  o2  «  SolEjiEj/Je.  The  only  apparent  difficulty  is  that,  unlike  o(e),  ofe^.-e)  may 

be  kinked  at  c  «  eJf  rendering  OjtEjjE)  discontinuous  at  this  point.  Examples  of  the 
two  exhaustive  alternatives i 

(2.19)  Alt.  1:  0(E  |E  )  =  o.  ,  Alt.  2:  o(E,;E_)  >  o,  , 

U  J  J  J  J  J 

are  shown  in  figure  3.  No  kink  occurs  for  Alt.  2;  indeed  in  this  case  o(c  ;•)  =  a ( * ) 
and  the  original  construction  (2.11)  -  (2.12)  yields  an  e  ( X )  obeying  (a)  -  (d).  The 


*or,  if  one  prefers,  o(£  ;e) 

J 


on  e  _  <  e  <  e 

J  u 
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Then  the  solution  <j>  8  S  of  (5.2)  ts  continuously  differentiable  with  $(0)  =  1  and 

O 

0  <  $(t)  <  1  for  all  t  >  0. 


The  theorem  will  be  proved  in  several  stages.  The  hypotheses  of  Theorem  2  are 
assumed  throughout. 

Lemma  1 :  Let  A:  B  +  R+ ,  where  the  set  B  C  R+  is  such  that  $(B)  C  B  and  0  e  B. 
Let  A  obey 

(5. 16)  0  <  A(t)  <  1 

and 


(5.  17) 

for  all  t  e  B. 

Proof :  If 

such  that  A(t 

kl 

then  Yg  =  1. 


[A(<p(t)>  -  2  -  Q(<Mt)  )  ]A(t)  +  1  =  0 
If  tk  e  B  and  t,  +  0,  then  Alt^)  ♦  1. 

A(  tu )  1  we  conclude  from  (5.16)  the  existence  of  a  subsequence  t. 

)  +  Yg  <  0  <  Yg  <  1.  Thus  it  is  sufficient  to  show  that  if  Att^)  +  'Ig' 


For  each  n  »  0,1,2,...,  (5.3)  yields 

(5.18)  lim  4>[nl(t  )  =  0  ,  lim  C(*lnl(tk>>  =  0  . 

k*“  k+°° 

Since  <p(B)  C  B,  (5.17)  ensures  the  existence  of  each  limit 

(5.19)  Y  *  lim  A(<^nj(tk>)  ,  n  =  1,2,..., 

k+* 

obeying 

(5.20)  Y_  “  2  -  1/Y  .  ,  n  =  1,2,..., 

n  ni 

which  is  the  same  as 

(5.21)  ,  -  =  (1  -  Yn_ ,  )/Y„_ ,  ,  n  =  1,2,... 

From  (5.19),  (5.16)  and  the  definition  of  Yq/  we  have  0  <  Yp  *  1  (n  =  0,1,...),  thus 

(5.21)  gives  Y  <  Y  .  .  Hence  y  ♦  Y  with  1  -  y  -  (1  -  Y)/Y,  so  Y  =  1.  Collecting 
n  n- 1  n 

results  we  have  1  =  Y  =  lim  y  <  Yq  *  1  /  whence  Yq  =  1« 

n*«> 
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for  all  n  =  1,2,3. . . 
(5.8)  we  arrive  at 
(5.12) 


hence  i|> 


and  all  t  >  0.  Taking  the  limit  of  (5.11)  as  n  +  »  and  using 

1 4»(t)  -  *<t>  I  <  | <t> ( o )  -  *(0)1  =  o  , 


o  o 

Theorem  1  implies  that  *  exists  a.e.  on  R+  with  0  <  $  4  1  and  moreover  that 

t. 

*(t)  =  /  *(s)ds,  the  integral  being  in  the  sense  of  Lebesgue.  Our  next  aim  is  to  show 
0 

that  *  is  continuously  differentiable  and  satisfies  (4.21).  Let  D  c  R+  be  a  Lebesgue 

°  c 

measurable  set  with  u(D)  =  0  such  that  *  exists  for  all  t  e  D  .  Here  and  throughout 
the  section,  u  will  denote  Lebesgue  measure  and  c  will  denote  set  complement  in  It*. 
Differentiating  (5.2)  gives 

o  o  o 

(5.13)  [*(*<t))  -  2  -  y(*(t))]*(t)  +1=0, 

for  all  t  e  Dc  such  that  *(t)  6  Dc.  Since  t  e  Dc  need  not  imply  *(t)  e  Dc,  it  is 
inconvenient  to  work  with  the  sets  D  and  Dc.  Instead,  we  seek  a  set  F  C  R+  such  that 
(i)  u(F)  =  0,  (ii)  *  exists  for  all  t  e  Fc  and  (iii)  4>(FC)  C  Fc.  clearly  the  set 


(5.14)  F  =  U  If*'1)  Inl  (D)}  , 

n=0 

Q 

satisfies  the  latter  two  requirements.  Thus  (5.13)  holds  for  all  t  e  F  .  To  show  that 
u(F)  *  0  we  note  that  )  8  S  ensures  that  ♦  ^  is  a  strictly  increasing  function. 
Moreover,  (5.2)  gives  explicitely 

(5.15)  4>~ 1  ( t )  =  2t  +  Q(t)  -  *(t>  . 

Suppose  Q  is  absolutely  continuous.  Then  so  is  *  ^  which,  in  turn,  renders 
u ( F )  =  0.  In  particular,  Q  is  absolutely  continuous  under  the  hypotheses  of 

Theorem  2:  Let  Q  be  continuously  differentiable  with  £(0)  =  Q(0)  =  0,  §(z)  >  0  for 


z  >  0  and,  in  addition,  let  5  be  Lipschitz  on  every  interval  [0,Zp]. 


He  are  first  going  to  show  that  the  iterative  scheme 

(5.7)  <t>  (t)  »  K-1(I4  ,(t)),  ♦,(t)  «  K-1(4  t)  , 

n  n- »  12 

converges  pointwise  to  a  solution  $  e  S  of  (5.5). 

We  have  0  <  $n(t)  <  since  $n  e  S.  An  Inductive  argument  now  gives 
#n<t)  >  ♦n-1(t),  thus  ♦  <M t )  for  all  t  >  0.  Since  $  e  S,  it  follows  that 

$(t)  is  (1)  Llpschitz  with  Lip[$]  <  1,  (ii)  increasing,  and  (iii)  obeys  iji(0)  “  0. 
Since  L$n<t)  ♦  I4(t)  and  K$n(t)  ♦  K$(t)  it  is  immediate  from  (5.7)  that  g>  is  a 
solution  of  (5.5)  or  equivalently  (5.2).  From  (5.2)  we  have  t1  *  t2  implies 
<Mt.,)  *  ♦(tj),  whence  $  is  strictly  increasing  and  so  $  e  S.  In  addition,  for 
t  >  0  we  have  pit)  >  0,  so  that  K-1(LiMt))  <  Li)>(t)  <  t  by  (5.4).  Since 
K  *(L$(t))  “  ♦(t)  ,  this  establishes  (5.3).  It  remains  to  show  that  $  is  the  unique 
solution  of  (5.2)  in  S. 

It  is  immediate  from  (5.3)  that 

(5.8)  lim  ♦lnl (t)  -  0  ,  (t  >  0)  . 

n-** 

Suppose  now  that  ip  e  S  is  also  a  solution  of  (5.2). 

Then 

I  <t>(  t)  -  t(t)|  -  |K-1(LiMt))  ~  K-1(LiKt))|  < 

(5.9) 

|  L$(t)  -  M>(t)|  -  j  |*(t(t))  -  ip(ip(t )  )  | 

<  i  |  <M<Ht) )  -  ip(4>(t ) )  |  +  j  | ip ( ( t  > )  -  >p<ip( t ) ) ) 

<  |  | ♦( ♦( t)  -  ip ( <P( t ) ) |  +  1  |*(t)  -  i|>(t)  |  , 

whence 

(5.10)  1  <P(t )  -  >p(t)  I  <  I <p ( <t> ( t > )  -  Ip ( 4> ( t )  )  I  . 

This  gives,  by  iteration, 

(5.11)  1 4>(  t)  -  p(t)|  <  l<p  (<p tnJ  (t>)  -  ip(<p [nI  ( t )  >  |  , 
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5.  Existence  of  a  Solution  to  the  Functional  Equation 


In  this  section,  we  prove  that  (4.14),  (4.19)  has  a  unique  solution  obeying  (4.21). 
In  addition,  the  first  of  the  two  theorems  in  this  section  furnishes  an  iterative  scheme 
for  solving  (4.19).  n 

Let  R+  *  {x|x  >  o}  and  for  n  =  0,1,2,  ...  let  $^(x)  =*  $($($(* >.($(x) )...)))  , 
♦  ^  (x)  =  x  .  Define  the  set  of  functions 

(5.1)  5  «  {♦sR*  ♦  R+|$(0)  *  0  ,  $(t)  strictly  increasing 

and  $  Lipschitz  with  Lip($]  <  l}  . 

Note  that  S  is  closed  under  composition.  The  following  theorem  shows  that  (4.19)  has  a 
unique  solution  $  e  S  for  a  class  of  functions  Q  which  include  those  Q  obeying 
(4. 14). 

Theorem  1 ;  Let  Q:  R+  ♦  *+  be  continuous  and  strictly  increasing  with  Q(0)  «  0.  Then 
the  equation 

(5.2)  ♦(♦It))  -  2*(t)  +  t  -  Q(<Mt>)  -  0  ,  (t  >  0) 
possesses  a  solution  $  e  5  .  Moreover  this  solution  <£  obeys 

(5.3)  0  <  <Mt)  <  t  ,  (t  >  0) 
and  is  the  only  solution  of  (5.2)  in  S. 

Proof :  Define  L:  S  +  S  by  Li(i(t)  »  <|>(i|i(t))  +  ~  t.  Also  define  K:  R  ♦  F 

by  K(z)  »  z  +  V2  Q(z) •  The  hypotheses  on  2  ensure  that  K  is  invertible  on  R+  with 

the  inverse  function  K-1  obeying 

(5.4)  0  <  K-1(z)  <  z  ,  <z  >  0)  . 

Furthermore,  one  may  verify  that  K-1  e  S. 

We  may  write  (5.2)  either  as 

(5.5)  K(  $  ( t )  )  -  L$(t)  ,  (t  >  0)  , 

or  as 

(5.6)  *(t)  -  K_1(L$(t>>  ,  (t  >  0)  . 
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Then  (4.15)  requires  that 

(4.17)  b($(z) )  =  -i  <J>  (  z )  +  jz  , 
whereas  (4.11)  demands  that 

(4.18)  b(4>(  z) )  =  G  ( <j>  (  z )  )  -  -1  <fr(4>(  z) )  -  -j  ♦(z),  Cb<  4>(Z) )  <  t)  . 

Consequently  (4.12),  (4.17-18)  yield  the  equation 

(4.19)  $(*<z))  -  2<j>(z)  +  z  =  Q(  4>  ( z )  )  , 

which  is  (1.5).  In  light  of  (4.16)  and  (4.9),  the  phase  boundary  trajectory  s(t)  is 
given  in  terms  of  $  as 

(4.20)  s(t)  -  c[a-1(t)  -  tj  ,  a(t)  =  1  [4(t)  +  tj  ,  t  <  t  , 

where  the  restriction  t  <  t  follows  from  b(4(z))  <  t.  Physically  for  t  >  t, 

E(s(t),t)  <  Cj  and  (4.8)  no  longer  reduces  to  (4.10).  Note  that  s(t)  completely 
determines  the  region  0  for  t  <  t  and,  by  virtue  of  (4.7),  the  displacement  u(x,t)  in 
this  region. 

The  derivation  of  (4.19)  assumes  that  s(t)  is  continuously  differentiable  and  obeys 
(4.5).  These  conditions  are  equivalent  to  $  being  continuously  differentiable  with 

O  O 

(4.21)  4(0)  -  0,  4(0)  =  1 ,  0  <  4 (z)  <1  for  z  >  0  . 

Notice  that  (4.21)  and  the  second  of  (4.20)  yields  a(t)  >  V2  »  which,  in  turn,  guarantees 
the  existence  of  a-1(t)  which  appears  in  the  first  of  (4.20). 


-19- 


h(t)  generates  acoustic  waves  which  drive  the  phase  boundary  forward  while  ringing  back 
and  forth  between  it  and  the  loading  device.  I  am  unaware  of  a  general  mathematical 
theory  appropriate  for  treating  this  system  in  the  above  form.  Fortunately,  this  problem 
simplifies  considerably  in  the  event  that  £(•»•)  and  v(*,*)  are  individually 
constant.  Recall  that  such  a  situation  prevails  immediately  ahead  of  the  phase  boundary 

»  A 

for  all  times  t  <  t.  Consequently,  for  all  times  t  such  that  b(t)  <  t,  one  has 
£(s(b(t) ) ,b(t) )  -  t(s(a(t) ),a(t) )  “  e  ,  v< s(b{ t) ) ,b(t) )  -  v(s(a(t) ),a(t) ) .  This,  in 

J 

conjunction  with  [o(£j)-dJ/c2  *  £j,  reduces  (4.8)  to 

(4.10)  h(t)  -  j  [£j  -  EjUSft)  -  a(t)  ]  +  £j,  (b(t)  <  t)  . 

This  equation  integrates  immediately  and,  with  the  aid  of  a(0)  -  b(0)  «  0,  may  be  cast 
into  the  form 

(4.11)  a(t)  +  b(t )  -  G(t) ,  (b(t)  <  t)  , 

where 

1  -  t  - 

(4.12)  G(t)  «  2t  ♦  jC(t),  Q(t)  -  [4/Uj  -  Ej))/  [h(r )  -  EjJdr  . 

We  note  for  future  reference  that  Q  is  of  unlimited  smoothness  with 

(4.13)  5(i)  -  4[h(z)  -IjJ/lEj  -  Ej)  , 
which  in  conjunction  with  (3.3)  yields 

(4.14)  Q(0)  «  5(0)  -  0,  5(z)  >0  (z  >  0)  . 

Return  now  to  (4.9)  and  note  that  (4.5)  ensures  that  a(t)  and  b(t)  are  each 
continuously  differentiable  and  monotonically  increasing.  Let  a  1(»),  b  ^ ( •  )  denote 
their  respective  inverse  functions.  Eliminating  s(»)  between  (4.9)  furnishes 

(4.15)  a-1(t)  +  b_1(t)  -  2t  . 

We  introduce 

(4.16)  ^(z)  »  2a(z)  -  z  ,  z  >  0  . 
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cue  may  express  u(s(t#),ta)  and  u(s(tb),tb)*  as  line  integrals  of  E(s(n)-»n), 
v(s(n)-,n)  along  £  »  s(n).  Upon  subsequently  eliminating  e(s(n)~,n),  v(s(n)-,n)  in 

favor  of  £(s(n)«h)»  v(s(n),n)  by  means  of  (1.4)  and  (4.1)  one  arrives  at 

u(x,t)  »  u( 0,0 )  +  /  *  {s(n)£(s(n) .n)  +  v(s(n).n))dn 

(4.7)  t  0 

-■5—  /  b  {[s(n)-c]v(s(n) .n)  -  cs<n)e(s(ii>  #n)  +  o(£(s(n) .n) )-D}dn 
t 

a 

The  function  s(n)  appearing  in  (4.7)  is  as  yet  unknown)  it  must  be  chosen  so  as  to 
satisfy  the  remaining  condition  e(0,t)  =  h(t).  With  a  view  toward  expressing 

e(0,t)  in  terms  of  quantities  along  £  «  s(n),  we  differentiate  (4.7)  with  respect  to 
x  and  subsequently  set  x  ■  0.  This  exercise  essentially  reduces  to  a  calculation  of 
3t  atfa 

— —  I  and  - —  I  .  These  may  be  expressed  in  terms  of  a(t)  and  £(t),  where 
3x  '  x**0  3x  1  x-0 

a(t)  i  ta(0,t)  and  b(t)  =  tb(U,t).  The  result  of  these  calculations  is 
h(t)  - -j  {[a(e(s(b(t) ) ,b(t) ) )-D  j/c2  -  e(s(b( t) > ,b(t) ) }b( t ) 

(4.8)  +|  {[a(?(s(a(t)),a(t)))-D]/c2  -  £(s(a(t) )  ,a(t> )  }a(t ) 

+  Y  |f(s(a(t)),a(t))  +  e  ( s(b(  t )  )  rb(  t ) )  +  [v(  s(  a(  t )  )  ,a(  t )  ) 

-  J(s(b(t)),b(t))]/c2}  . 

The  functions  a<t)  and  b(t)  are  connected  to  s(t)  by  means  of  the  implicit  relations 

(4.9)  s(a( t) )  -  c(t  -  a(t)]  -  0,  s(b(t))  +  c[t  -  b( t ) ]  =  0,  a(0)  -  b(0)  -  0  . 
Equations  (4.8),  (4.9)  for  the  unknowns  a(*)»  a(*),  b(*)  describe  how  the  end-strain 


‘which  are  continuous  across  £  -  s(n)  by  virtue  of  the  first  of  (1.4). 
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values  behind  the  phase  boundary  by  means  of  the  discontinuity  conditions  (1.4).  With  the 


aid  of  (4.1),  these  discontinuity  relations  can  be  solved  for 

o 

v(s(n)“»h)  provided  that  s(n)  *  c.  In  the  interior  of  U, 
(1.1),  which  becomes 


(4.4) 


2 

c 


e(s(n>-,n)  and 

the  dynamical  fields  obey 


by  virtue  of  (4.1). 

We  assume,  subject  to  later  verification,  that  s(n)  is  continuously  differentiable 


with 

(4.5)  s(0)  *=  s(0)  =  0  and  0  <  s(n)  <  c  for  n  >  0  . 

The  standard  representation  for  solutions  of  (4.4)  in  terms  of  Cauchy  data  on  a  given 
curve  furnishes 


(4.6) 


u(x,t) 


»  V2  {u(s(ta),ta)  +  u(s(tb),tb)} 

-■5^  /  b  {v(s(n)-,n)s(n)  +  c2E(s(n)-,n) }dn 

fca 


for  all  ( x, t )  e  0.  Here,  ta  «  ta(x,t)  and  tjj  *  t^(x,t)  delimit  the  domain  of 
dependence  upon  £  «  s(n)  as  depicted  in  figure  4. 


Figure  4.  Graphical  representation  of  the  times  tb  =  t^(x,t)  and  ta  =  ta(x,t) 

for  (x,t)  e  U.  These  times  are  also  given  implicitely  as  the  roots  of 
x  »  s(tb)  ♦  c[t  -  t^)  and  x  =  s(ta)  -  c[t  -  ta)  respectively. 


4.  A  Functional  Equation 


For  the  remainder  of  this  paper  we  restrict  attention  to  materials  for  which 
i)  aa  <  og, 

(4.1) 

2  — 
ii )  o(e)  =  c  e  +  D  ,  e  >  tj  . 

Such  "second  phase  linear"  materials  have  been  studied  previously  in  [11]  in  connection 
with  this  problem  in  which  Oj  =  og.  We  shall  show  that  there  exists  a  time  t  >  t*,  and 
possibly  infinite,  such  that  P(t)  has  a  solution  with  an  energy  confining  phase 

boundary,  the  location  of  which  is  governed  by  (1.5). 

* 

By  transforming  the  time  variable  t  +  t  -  t  we  may,  without  loss  of  generality, 
set  t*  »  0.  Let  (£,n)  be  a  point  in  the  (x,t)-plane  and  let  V  and  U  be 
respectively  the  regions  in  front  and  behind  the  phase  boundary.  In  terms  of  (£,n),  we 


v  =  {(5,n)  |  5  >  s(n),  n  >  o|  , 

(4.2) 

v  -  {(£,n)  |  o  <  5  <  s(n),  n  >  o}  . 

The  fundamental  difference  between  U  and  V  is  that  whereas  the  dynamical  fields 
e  and  v  in  0  remain  to  be  found,  these  fields  in  V  are  determined  by  events 
occurring  at  x  =  0  before  the  time  of  phase  boundary  emergence.  To  acknowledge  this 
distinction,  we  shall  let  e  and  v  be  respectively  the  (known)  strain  and  velocity  in 
the  region  V.  In  particular,  £  and  v  are  individually  constant  in  that  portion  of 
V  between  the  phase  boundary  5  =  s(n)  and  the  characteristic  curve  £  c  /o' ( Ej)  n* 
Should  the  phase  boundary  overtake  the  front  of  this  region,  we  shall  denote  by  t  the 
time  at  which  this  event  occurs.  Thus  t  >  0  is  the  (first  positive)  root  —  if  it 
exists  —  of 


(4.3)  s(t)  =  /o' (e  >  t  ,  t  >  0  . 

J 

We  shall  say  that  t  =  “  in  the  event  that  the  phase  boundary  never  overtakes  this  front. 

For  all  values  n  >  0  the  strain  and  velocity  immediately  in  front  of  the  phase 
boundary  —  e(s(n)+,n)  =  e(s(n),n)  and  v(s(n)+,n)  =  v(s(n),n)  —  are  functions  of  the 
values  n  and  s(n).  This  strain  and  velocity  is  in  turn  related  to  the  corresponding 


particular,  tha  characteristic  C+  which  would  issue  from  (x,t)  -  (0,t*)  is  the  ray 


x  -  /o'Uj)  [t  -  t*l. 

At  the  time  t  «  t*  a  phase  boundary  —  x  -  s(t),  say  —  emerges  at  x  «  0  and 
subsequently  propagates  into  tha  interior.  By  virue  of  (1.4),  the  phase  boundary  velocity 
obeys 


(3.4) 


s(t) 


r  q(e(s(t)+,t) )  -  o(e(»(t)-,t) )  i  V2 
'  E(s(t)+,t)  -  e(s(t)-,t)  • 


This  relation,  in  conjunction  with  the  continuity  of  g(t)  at  t  «  t*,  yields 
s(t*)  -  0.  Thus  there  exists  a  region  in  the  (x,t)~  plane 

(3.5)  s(t)  <  x  <  /o'(e,)  [t  -  t*]  , 

immediately  in  front  of  the  phase  boundary,  which  is  conspicuous  because  the  solution 
within  this  region  is  not  as  yet  determined.  We  now  inquire  as  to  the  value  of  the  strain 
in  this  region. 

A  first  alternative  is  that  the  value  of  £(x,t)  in  this  region  could  exceed  the 
value  Cj  and  yet  remain  below  the  value  81 •  If  such  is  the  case,  then  events  at  the 
external  boundary  occurring  after  t  “  t*  should  play  some  role  in  determining  £(x,t) 
in  this  region.  A  second  alternative  —  that  which  shall  be  pursued  here  ~  is  that  the 
strain  in  this  region  remains  fixed  at  the  value  £j.  In  this  case,  the  velocity  in  this 
region  must  also  be  constant,  its  value  being  given  by  the  (constant)  velocity  on  the 
boundary  characteristic  x  «■  /o'(e,)  (t  -  t*J.  Uhlike  the  first  alternative,  this 
situation  is  one  in  which  events  at  the  external  boundary  occurring  before  t  «  t* 
completely  determine  the  fields  in  front  of  the  phase  boundary.  What  is  not  determined  is 
the  location  of  this  boundary  itself,  as  well  as  the  fields  behind  it.  Physically  this 
corresponds  to  a  situation  in  which  all  of  the  energy  delivered  to  the  bar  after  the 
emergence  of  the  phase  boundary  remains  confined  behind  the  phase  boundary.  We  shall  use 
the  expression  energy  confining  phase  boundary  for  this  second  state  of  affairs.  In  the 
next  section  we  formulate  a  free  boundary  problem  for  such  a  phase  boundary  and  show  that 
it  leads  to  equation  (1.5). 


3.  Energy  Confining  those  Boundaries  in  the  Presence  of  Smooth  Loads 


Me  now  focus  attention  upon  loads  g(t)  which  are  of  unlimited  smoothness  with 
g(0)  »  0,  g(t)  >  0  for  0  <  t  <  t*  and  g(t)  >  Oj  for  t  >  t*.  It  will  be  convenient 
to  write  the  boundary  condition  (1.2)  ae 

(3.1)  e(0,t)  -  h(t)  » 

where  h(t)  is  the  end-strain,  obeying  (2.1),  associated  with  the  traction  g(t).  It  is 
given  uniquely  as  the  solution  to 

(3.2)  g(t)  -  oUjihtt))  . 

The  function  h(t)  of  (3.2)  is  discontinuous  at  t  -  t*,  while  at  all  other  values  t  it 
is  of  unlimited  smoothness.  Specifically  it  obeys 

(3.3)  h( 0)  -  0,  fi(t)  >0  (0  <  t  <  t*), 

h(t*-)  -  e .  h(t*+)  -  7_,  h(t)  >  (t  >  t*)  . 

J  J  J 


say  £(£  ,X)  where  e(e(X^+),X)  »  c(X).  All  the  solutions  for  £*  <  e ( X ^  )  have  phase 
boundaries  which  travel  faster  than  that  of  the  solution  £(X). 

Now  suppose  that  e(X^-)  <  e Then  the  upper  concave  envelope  of  0(e)  deleted  on 

(£J,£**)  with  e ( X ^— )  <  £**  <  again  differs  from  office)  and  the  resulting 
solutions  all  have  slower  moving  phase  boundaries  than  that  of  the  solution  e(X). 
Solutions  with  more  than  N  discontinuities  can  result  from  deleting  o(e)  on  additional 
intervals*  Finally,  the  extreme  case  of  deleting  everything  but  the  end-points  e  »  0 
and  e  «  yields  an  upper  concave  envelope  which  is  the  line  segment  connecting 
(0,0)  to  l  c„*  <>„,)•  3*>e  associated  solution  is 


•  0  » 
£  (  X)  - 


(2.26) 


x  >  . 

0  <  X  <  /o  /£  ' 

00  00 


X  >  ^o«/e-  « 

o  <  x  <  /7jT„  . 


As  is  the  case  with  £,  these  other  solutions  e  do  not  satisfy  all  of  the  inequalities 

(2.13)  -  (2.16)  unless  f(e)  -  0(e)  for  all  0  <  e  <  e^. 

In  summary,  whenever  o(e  >e)  differs  from  the  line  segment  connecting  (0,0)  to 

J 

(e^.a^l,  there  exist  solutions  other  than  e(X),  v(X)  for  the  impulsive  load  (2.3).  In 

particular,  a  family  of  solutions  occurs  for  Alt.  1,  since  then  e(X+)=E>0.  A 

1  J 

pertinent  question  is:  which,  if  any,  of  these  solutions,  gives  an  asymptotic  description 
of  the  dynamical  fields  in  a  bar  subject  to  a  smooth  load  g(t)  with  g(t)  oj  In  what 
follows  it  will  be  shown  under  certain  conditions  for  a  special  class  of  materials  that 


the  problem  with  a  smooth  load  is  reducible  to  equation  (1.5).  One  result  of  examining 
this  equation  asymptotically  will  be  that  if  g(t)  ♦  and  Alt.  1  holds,  then  £(x/t) 
is  the  x/t  fixed,  t  +  »  limit  of  the  resulting  strain  field. 


n 


Finally,  e(X)  obeys  (2.5)  whenever  .»  graph  of  o  (e  ;e)  is  neither  horizontal  nor 
vertical.  Thus  extending  e(X)  to  X  >  o  by 


(2.22) 


c  (  X )  =  0  , 


e(X)  =  £.  , 


X  >  Aj_(e  >0)  , 

Z  J 


0  <  X  <  So2(eJi£xi)  , 


yields  an  c ( X )  which  fulfills  all  the  requirements  (a)  -  (d). 

An  interesting  feature  of  this  solution  is  that  the  phase  boundary  may  travel  slower 
than  the  sound  speed  of  the  material  into  which  it  is  proceeding.  Although  (2.13),  (2.15) 
remain  in  force  for  ordinary  shocks,  (2.14),  (2.16)  continue  to  hold  only  for  Alt.  2.  For 
Alt.  1,  (2.14)  is  replaced  by 


(2.23) 


»'(e(X  -))  -  X,  <  /oMElX^))  -  /o'U  )  ,  (N  >  1,  Alt.  1) 


whereas  (2.16)  is  changed  to 


(2.24) 


X^  j,  <  /o'(e  ( X1 +) )  ■  /o'IEj)  ,  (N  -1,  Alt.  1) 


We  now  reserve  the  notation  €(X)  for  that  particular  solution  of  (a)  -  (d)  given  by 
(2.17),  (2.18),  (2.22)  and  enquire  into  the  uniqueness  of  this  solution. 

The  construction  leading  to  e(X)  succeeded  in  satisfying  (d)  while  maintaining 
(a)  -  (c)  by  "deleting"  0(e)  on  <  t  <  e^  befo-e  taking  the  upper  concave 

envelope.  Indeed,  if  even  more  of  0(e)  is  deleted,  giving  say  f(e),  then  the  ensuing 

*  m 

concave  envelope  f(e)  leads  to  a  solution  e(A)  of  (a)  -  (d)  by  the  same  procedure 


applied  to 


(2.25) 


X  »  f  ,  e  (  X )  j 


Thus,  for  example,  if  €  ( X  ^  +- )  >  0,  then  the  upper  concave  envelope  of  alt)  deleted  on 

(e*,7, ),  with  0  4  e*  <  e(A  +),  differs  from  o(e  ;£)  and  the  resulting  e(X)  differs 
J  I  J 

from  £ ( X ) .  Moreover  each  choice  of  £*  on  the  interval  0  <  £*  <  e(X^+)  generates  a 
different  upper  concave  envelope  t(A)  so  that  t*  parameterizes  a  family  of  solutions. 


-  1 1- 


c(e  .*£) 


Figure  3.  The  upper  concave  envelope  o{e_/e)  for  two  different  values  of  e  . 

J  00 

This  envelope  need  only  be  above  the  dashed  curve,  which  is  o(e_;e). 

J 

The  left  and  right  figures  are  examples  of  Alt.  1  and  Alt.  2  respectively, 
program,  however  must  be  modified  for  Alt.  1,  because  it  is  generic  for  this  alternative 


(2.20) 


o'Uj)  *  >  °2lEJ,ej+> 


In  the  face  of  (2.20)  let  us  agree  to  identify  °2*ej,ej*  with  the  full  interval 

[°2  ^  ^ '  ^2  ^  graph  of  will  then  be  continuously  decreasin 

on  0  <  e  <  and  will  involve  a  vertical  segment  at  e  =  eT.  Hence  (2.18)  defines  a 

J 

piecewise  continuous  e(A)  on  Jo  At  ;c  )  <  A  <  Jo ,  (eT;0>  .  if  A,  >  A.  >  ...»  A 

i  j  12  N 

are  the  values  where  £(•)  is  discontinuous,  then  A^  is  the  height  of  the  i  — 

A  A 

horizontal  portion  of  o  (e  >e).  On  the  other  hand,  the  vertical  segment  of  o,(e  >e )  a 

*  d  2  J 

e  *  e_  yields 

d 


e(A)  *  e,  ,  /  a. (e , ; e,+)  <  A  <  /a ' ( £ _ > 


r\ 


Two  sets  which  satisfy  the  requirements  upon  b  stated  above  are  F  and  R  •  Note 

o 

that  $  obeys  (5.16),  (5.17)  for  B  =  Fc.  Thus  we  draw  three  corollaries  central  to  the 
following  development 

Corollary  1 ;  All  A:  FC  ♦  R+  obeying  (5.16),  (5.17)  on  Fc  have  the  property 

that  A( t  )  +  1  whenever  t,  ■*  0,  t,  e  Fc. 
k  k  k 

O 

Corollary  2:  $(t^)  +  '  for  +  0,  t^  e  Fc. 

Corollary  3;  If  As  R+  *  R+  obeys  (5.16),  (5.17)  on  R+,  then  A  is 
continuous  at  the  origin  with  A(0)  *  1. 

O 

The  first  two  corollaries  guarantee  that  ^  is  the  unique  solution  of  (5.17)  in  the 
following  sense 

Lemma  2 :  If  A:  R+  ♦  R+  obeys  (5.16),  (5.17)  for  all  t  e  Fc,  then 
Aft)  *  $(t)  for  all  t  e  Fc, 

Proof :  write  (5.17)  in  the  form  A(t)  =  -1/[A($(t))  -  2  -  &($(t))];  a  corresponding 

O 

result  holds  with  $  replacing  A.  By  subtracting  these  two  expressions  and  invoking 
O  o 

(5.16),  0  k  $  <  1 ,  and  Q(<|>(t))  >  0,  we  arrive  at 


(5.22) 


I  A( t)  -  <fr(t)|  k  ] A( $( t ) }  -  *(*(t))|  ,  (t  e  r  )  , 


c  c 

which  gives,  by  induction  since  ♦(F  )  C  F  , 


(5.23) 


|  A ( t )  -  <f>(  t)  |  <  |  A(  <|) [n]  <  t )  ) 


(t))|  ,  (t  e  FC,  n  »  1.2,...) 


Letting  n  we  have  from  (5.8)  and  Corollaries  1  and  2, 

(5.24)  |  A(  t)  -  <(>(  t )  |  k  lim  |  A(  r )  -  *°(r)|  «  0  ,  (t  e  Fc)  ./ 

r+0.  f 


We  already  have  that  $  is  a  solution  of  (5.17)  on  Fc.  We  now  show  that  (5.17) 
possesses  a  solution  not  only  on  Fc,  but  also  on  R+. 

Lemma  3;  There  exists  A:  R+  ♦  R+,  continuous  with  A(0)  *  1,  0  <  A(t)  <  1 

for  t  >  0  such  that  (5.17)  holds  for  all  t  e  R+. 


Proof ;  Consider  the  sequence  of  functions  generated  by 


( 5.25)  \+1(t)  “  V{2  +  QUU))  -  A^tt))}  ,  A^t)  «  1/{  2  +  ii(*(t))}  . 

O 

Separate  arguments  based  upon  induction  and  Q(z)  >  0  yield  0  <  A^(t)  <1,  (t  >  0), 

and  A^+1<t)  >  A^ft),  (t  >  0).  Hence  Aft)  =  Xim  A.  (t)  exists  and  obeys 

k+*» 

0  <  Aft)  <  1.  It  is  evident  from  (5.25)  that  Aft)  is  a  solution  of  (5.17).  Also 

r 

A^l)  «  1/{2  +  Q(4>( 0) )  -  A( $( 0) )  }  -  1/{2  -  A(0 ) }  so  that  A(0)  -  1.  Moreover,  for 

O  O 

t  >  0,  the  strict  inequality  Qf^ft))  >  0  furnishes  1/A(t)  -  2  +  Q ( <#> ( t ) )  -  A< <f>  ( t  > )  > 

2  -  A($( t) )  >  1,  whence  0  <  Aft)  <  1  for  t  >  0.  The  continuity  of  A  at  the  origin 
is  the  result  of  Corollary  3.  It  remains  to  show  that  Aft)  is  continuous  for  t  >  0. 

Let  t  >  0  and  e  >  0  be  given)  we  shall  show  that  there  exists  6  >  0  such  that 
if  r  >  0  with  | t  —  t |  <5  then  |A(t)  -  A(t)|  <  e.  Let  q  be  a  Lipschitz  constant 

O 

for  Q  on  0  <  z  <  t  +  1.  Choose  d  >  0  such  that  0  <  x  <  d,  0  <  y  <  d  implies 

|a(x)  -  A ( y ) j  <  e/2.  Now  (5.8)  guarantees  the  existence  of  an  N  such  that 

INI  ° 

J  ft  +  1)  <  d.  From  (5.17)  and  the  bounds  Afz)  <  1,  Qfz)  >  0,  it  follows  that 

(5.26)  |  Aft )  -  Aft)  I  <  Iflf  $(T) )  -  Q($<t))|  +  |A(*(t))  -  A(  <J.(  T  >  >  |  ,  ft  >  0,  T  >  0)  . 

Iterating  this  inequality  N  -  1  times,  we  arrive  at 

(5.27)  |  A(  t)  -  A(T)|  <  {  l  |Q<  tK1  ft))  -  Q(*[kl(T)>)}  +  |  A(  lN1  ft))  -  A(4>lN)  (T))|  . 

k-1 

Let  6  «  min{l,  e/2Nq}  and  suppose  that  | t  —  T |  <6.  Then  (5.27)  implies 


(5.28) 


|A(t)  -  A(T)|  <  {  l  q|t  -  t| }  +  |  <  e  . / 
k-1  ' 


we  now  rapidly  complete  the 


Proof  of  Theorem  2:  let  A  be  as  given  in  Lemma  3.  Then  Lemmas  2  and  3  furnish, 
since  u(F)  -  0, 


-25- 


6.  Large  Time  Behavior  of  Solutions 

Theorems  1  and  2  guarantee  the  existence  of  a  solution  $(t)  to  the  system  (4.14), 
(4.19),  (4.21).  The  function  in  turn,  enables  one  to  determine  the  region  U  and 

the  function  u(x,t)  in  (J  by  means  of  (4.6),  (4.20).  Thus  we  arrive  at  a  solution 
involving  an  energy  confining  phase  boundary  to  T(t)  for  the  class  of  materials  (4.1). 

A 

Here  t  is  found  from  (4.3).  Thus  the  question  arises,  what  circumstances  —  if  any  — * 
give  rise  to  a  problem  for  which  t  =  ®?  Or,  equivalently,  in  terms  of  the  known  fields 
in  front  of  the  phase  boundary:  what  conditons  give  rise  to  a  phase  boundary  which  never 
overtakes  the  front  of  the  constant-strain  region? 

Since  the  leading  edge  of  this  front  travels  at  speed  /o' (£j) ,  a  partial  answer  to 

the  above  question  is  immediate  from  the  bound  s(t)  <  c.  Namely,  if  /o' (£j)  >  c, 

then  s ( t )  <  Jo' ( e,)  which  yields  t  =  “.  On  the  other  hand,  for  materials  with 
J 

_  A 

Jo' ( e,)  <  c,  we  will  now  show  that  t  <  “  provided  h(t)  becomes  sufficiently  large. 

J 

Consider  h(t)  *  £  >  eT  as  t  ♦  Then  from  (4.12)  it  follows  that 

00  J 


(6.1) 


whence  (4.19)  yields,  to  leading  order  as  t  + 

(6.2)  $(<Mt))  +  t  ~  (2  +  4 

where  we  have  introduced  y  =  4[eoo  -  ej]/[e„  - 
[y  +2]$(t)  +  t  «  0  has  two  solutions,  ^(t)  ; 


■»,  the  asymptotic  relation 


}$(t)  =  [y  +  2  ]  *  (t )  , 

e,l  >  0.  The  equation  $($(t)) 
vJ  J 

J  [y  +  2  ±  Jy2  +  4y ]t .  However 


since 


Y  >  0,  only  the  negative  root  yields  a  $(t)  obeying  0  <  iji(t)  <  t  for  t  >  0.  Hence 


Assume  now,  for  the  sake  of  argument,  that  t  -  <*.  Then  (4.20)  leads  to 


>(t) 


- P--)  *-/S 

+  4  -  /y  +  4Y  / 


ct  , 


(6.4) 


s(t)  ~  j— 


e„  -  et 
"  j 


e»  *  ej 


In  view  of  (4.3),  (6.4)  a  necessary  condition  for  the  assumption  t  «  •  to  hold  is 


(6.5) 


<  /o' (Ej)  . 


This  condition  is  satisfied 


(6.6) 


for  e,  <  e„<  e„  if  /o'  (e  .  )  <  c  , 

J  ®  1  1 


for  all  >  tj  if  /o' <  )  >  c  , 


where  the  value 

(6.7) 


£jg,(ej) 

(£j) 


The  second  of  (6.6)  was  anticipated  from  the  previous  discussion.  The  first  of  (6.6) 

provides  a  boundary  data  upper  threshold  which  is  necessary  for  the  phase  boundary  to 

remain  behind  the  front  of  the  constant  Btrain  region;  hence  it  is  a  necessary  condition 

for  t  ■  We  do  not  enquire  in  detail  as  to  conditions  sufficient  to  ensure  that 

t  «  “  for  the  case  /o' ( e .)  <  c,  other  than  to  remark  that  it  is  evident  from  (4.3) 

J 

•  /'  *  1  A 

that  if  s(t)  <  /o'(e,)  for  all  t  >  0,  then  t  ■  “.  Thus  c  <  e„  is  also  a 
J  ®  T 

*  oo 

sufficient  condition  for  t  *  •  whenever  s(t)  >0  for  all  t  >  0.  We  conjecture  that  a 


condition  sufficient  to  preclude  a  decellerating  phase  boundary  is 


For  the  cases  h(t)  ♦  £„  >  £j  with  t  “  »,  we  may  calculate  the  large  time  strain  and 
velocity  fields  for  the  problem  P(“),  (4.1),  in  the  limit  t  *  »  with  X  ■  x/t  fixed. 

As  mentioned  in  Section  3,  these  fields  in  the  region  x  >  /o' (£j)  tt  -  t*] ,  are  calculated 
in  reference  [11].  Whence  we  verify  for  x  >  {[e^  -  ejJ/[e®  “  Ej]}  ct  ~  s(t),  that  e 
and  v  admit  the  expansions  £(x,t)  “  £(Xt,t)  «  e^X)  +  £2<X,t),  v(x,t)  - 
v^X)  +  v^tX.t)  where  e2,  v2  ♦  0  as  t  ♦  ■  for  fixed  X,  and 

(i)  e^X)  -  0,  v , (X)  -  0, 

(ii)  e^X)  -  [O*]”1  (X2)  , 
e^X) 

(6.8)  v.(X)  «  -  /  /o* ( s)  ds  , 

0 

( iii )  £ 1 ( X )  -  £j  , 

V  (X)  -  -  /  J  /oMbT  dl  :  V  , 

0 

Behind  the  phase  boundary,  (4.7)  leads  to  the  exact  expressions: 

£(x.t)  »  [<£j  -  v*/c)s(tb)  +  (“dj/c  +  v  +  D/c)]  /  [2(s(tfc)  -  c)] 

+  [ ( £ j  +  v*/c)s(ta)  +  (°j/c  +  v  -  D/c)]  /  [2(s(ta)  +  c ) ]  , 

(6.9)  v( x, t )  -  [ ( £ j  -  v*/c)s(tb)  +  (Oj/c  +  V  -  D/c)]  /  [2(1  +  s(tb)/c)] 

♦  [(Ej  ♦  V*/c)5(ta)  +  (Oj/C  +  V  -  D/C)]  /  [2(1  +  8(ta)/c)J  , 

where  v*  is  the  same  value  defined  in  (6.8)  and  ta  *  ta(x,t),  tb  «  tb(x,t)  are,  as 

before,  depicted  in  Figure  4.  Now  it  follows,  either  from  Figure  4  —  or  by  writing 
expressions  for  ta  and  tb  —  that  t#  ♦  •,  tb  ♦  “  in  the  limit  t  ♦  *°,  x/t  fixed, 

•  •  ”  J 

whence  s(ta)  ♦  s(tb)  ♦  [(e^  “  ej)/(eM  “  Ej)l  ‘c  Xn  this  same  limit.  Entering  (6.9)  with 
this  result  yields,  as  t  ♦  ",  X  *  x/t  fixed, 


/o' (0)  4  X  4  «  , 


/o*(E j)  <  X  <  /a’(O)  , 


li  s  / T~r V  <=  <  x  <  /o- 


(ej> 
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e(x,t)  ♦  tm  , 

(6.10) 

) - “ 

v(x,t)  ♦  -  /(c.  -  -  £j)  c  +  V  ,  0  <  A  <  A,  . 

For  materials  obeying  (4.1)  and  values  of  ea  subject  to  (6.6),  the  right-hand  sides 
of  (6.8),  (6.10)  furnish  an  exact  solution  of  the  form  (2.4)  to  the  problem  P (°“ )  for  the 
impulsive  load 

^0,  t  -  0  , 

(6.11)  h(t)  * 

's‘  >  eJ(  t  >  0  . 

Since  this  solution  also  satisfies  conditions  (a)  -  (d)  of  section  2,  we  shall  conpare 
this  solution  to  those  impulsive  load  solutions  obtained  in  that  section. 

For  the  materials  now  under  consideration,  condition  (6.6)  implies  (and  is  implied 

by)  Alt.  1  of  (2.19).  Recall  that  for  this  alternative,  the  impulsive  load  problem  has  a 

«  *  * 

one  parameter  family  of  solutions  e(e,A),  0  <  e  <  However,  of  these  solutions, 

only  e(e  ,A)  =■  e(A)  given  by  (2.17),  (2.18),  (2.22)  has  a  region  in  front  of  the  phase 
boundary  with  E  «  E  •  Indeed,  it  is  easily  verified  that  this  £(X)  is  given  precisely 
by  the  right  hand  sides  of  (6.8),  (6.10).  Thus  6(A)  given  by  (2.17),  (2.18),  (2.22) 
is  —  at  least  for  materials  obeying  (4.1)  and  loads  obeying  (6.6)  ■—  appropriate  for 
energy  confining  phase  boundaries. 

We  close  this  paper  on  the  remark  concerning  the  case  of  a  smooth  h(t)  *  >  eT 

for  materials  (4.1)  with  /o'(e,)  <  c.  Then,  as  shown,  t  <  *•,  and  the  solution  $  of 

J 

(4.19)  does  not  lead  to  a  solution  of  P (T)  for  T  >  t.  Nevertheless,  we  conjectue  that 
P(«*)  has  a  solution  involving  an  energy  confining  phase  boundary  and,  moreover, 
anticipate  that  E(A)  given  by  (2.17),  (2.18),  (2.22)  is  the  A  **  x/t  fixed,  t  +  • 

limit  of  the  strain  field.  In  particular,  this  predicts,  in  place  of  (6.4),  that 

1  /  » 

lim  s(t)  <  { [ -  E^]  /  [eb  -  E^j}  /2c.  The  problem  P(T),  T  >  t  may  be  studied  in  the 

framework  developed  here  by  treating  the  system  (4.8),  (4.9)  which  describes  an  energy 
confining  phase  boundary  advancing  into  a  more  general  strain  and  velocity  field. 


-30- 


Acknowledgement : 


The  author  would  like  to  acknowledge  helpful  discussions  with 
James,  w.  Hrusa,  R.  Lyons,  M.  Renardy  and  R.E.L.  Turner. 


-31- 


REFERENCES 


1.  D.  E.  Grady,  R.  E.  Hollenbach  and  K.  w.  Schuler,  Compression  Wave  Studies  in  Calcite 
Rock,  J  Geophys.  Res.,  83  (1978),  p.  2839. 

2.  R.  D.  James,  Co-Existant  Phases  in  the  One- Dimensional  Static  Theory  of  Elastic  Bars, 
Arch.  Rational  Mech.  Anal.,  72  (1979),  p.  99. 

3.  R.  D.  James,  The  Propagation  of  Phase  Boundaries  in  Elastic  Bars,  Arch.  Rational 
Mech.  Anal.,  73  (1980),  p.  125. 

4.  J.  R.  Yaeger,  A  Practical  Shape-Memory  Electromechanical  Activator,  Mechanical 
Engineering,  7  (1984),  p.  51. 

5.  J.  L.  Ericksen,  Equilibrium  of  Bars,  J.  Elasticity,  5  (1975)  3-4,  p.  191. 

6.  J.  K.  Knowles  and  E.  Sternberg,  On  the  Failure  of  Ellipticity  and  the  Emergence  of 
Discontinuous  Deformation  Gradients  in  Plane  Finite  Elastostatics,  J.  Elasticity,  8 
(1978),  p.  329. 

7.  J.  K.  Knowles  and  E.  Sternberg,  Discontinuous  Deformation  Gradients  Near  the  Tip  of  a 
Crack  in  Finite  Anti-Plane  Shear:  An  Example,  J.  Elasticity,  10  (1980),  p.  81. 

8.  J.  K.  Knowles  and  E.  Sternberg,  Anti-plane  Shear  Fields  With  Discontinuous 
Deformations  Gradients  Near  the  Tip  of  a  Crack  in  Finite  Elastostatics,  J. 

Elasticity,  11  (1981)  p.  129. 

9.  R.  c.  Abeyaratne,  Discontinuous  Deformation  Gradients  in  the  Finite  Twisting  of  an 
Incompressible  Elastic  Tube,  J.  Elasticity,  11  (1981)  p.  43. 

10.  R.  c.  Abeyaratne,  Discontinuous  Deformation  Gradients  Away  From  the  Tip  of  a  Crack  in 
Anti-plane  Shear,  J.  Elasticity,  11  (1981),  p.  373. 

11.  T.  J.  Pence,  On  the  Bnergene  and  Propagation  of  a  Phase  Boundary  in  an  Elastic  Bar 
With  a  Suddenly  Applied  End  Load,  to  appear  J.  Elasticity. 

12.  J.  Serrin,  Phase  Transitions  and  Interfacial  Layers  for  Van  der  Waals  Fluids,  in 
"Proceedings  of  SAFA  IV  Conference,  Recent  Methods  in  Nonlinear  Analysis  and 
Applications,  Naples,  March  21-28,  1980"  (A  Canfora,  S.  Rionero,  C.  Sbordone,  C. 
Trombetti,  Eds.) 


-32- 


13.  M.  Slemrod.  Admissibility  Criteria  for  Propagating  Phase  Boundaries  in  a  Van  der 
Waals  Fluid,  Arch.  Rational  Mech.  Anal.,  81  (1983),  p.  301. 

14.  R.  Hagan  and  M.  Slemrod,  The  Viscosity  -  Capillarity  Criterion  for  Shocks  and  Phase 
Transitions,  Arch.  Rational  Mech.  Anal.,  83  (1983)  p.  333. 

15.  M.  Kuczma,  Functional  Equations  in  a  Single  Variable,  PWN-Polish  Scientific 
Publishers,  Warsaw,  (1968). 

16.  M.  P.  Mortell  and  B.  R.  Seymour,  Exact  Solutions  of  a  Functional  Equation  Arising  in 
Nonlinear  Wave  Propagation,  SIAM  J.  Appl.  Math.,  Vol  30,  No.  4  (1976),  p.  587. 

17.  M.  E.  Gurtin,  On  a  Theory  of  Phase  Transitions  With  Interfacial  Energy,  preprint. 


TP/jp 


-33- 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Phan  Data  Bntarad) 


REPORT  DOCUMENTATION  PAGE 


REPORT  HUM 

#2827 


«.  TITLE  (and  Subtltla) 

A  FUNCTIONAL  EQUATION  GOVERNING  MOVING 
PHASE  BOUNDARIES  IN  AN  ELASTIC  BAR 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


S.  RECIPIENT’S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Summary  Report  -  no  specific 
reporting  period 


6.  PERFORMING  OR'  REPORT  NUMBER 


7.  AUTHORS 


B.  CONTRACT  OR  GRANT  NUMBERf*) 


Thomas  J .  Pence 


•  ■  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Mathematics  Research  Center,  University  of 
610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53706 


H-  CONTROLLING  OFFICE  NAME  AND  AODRESS 


I 


(See  Item  18  below) 


MONITORING  AGENCY  NAME  ft  ADDRESSfll  dlllarant  from  Controlling  Olll  co) 


IS-  DISTRIBUTION  STATEMENT  (ol  thla  Raport) 

Approved  for  public  release;  distribution  unlimited. 


DAAG29-80-C-0041 

MCS-8210950 


10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

Work  Unit  Number  2  - 
Physical  Mathematics 


12.  REPORT  DATE 

June  1985 


IS.  NUMBER  OF  PAGES 

33 


IS.  SECURITY  CLASS,  (of  thla  raport) 

UNCLASSIFIED 


mm 


'7.  DISTRIBUTION  STATEMENT  (of  111*  abatraet  antarad  In  Block  20,  If  dlllarant  from  Raport) 

phase  transitions,  elastic  solids.  Functional  equations 


IB.  SUPPLEMENTARY  NOTES 

U.  S.  Army  Research  Office  Natii 

P.  0.  Box  12211  Wash; 

Research  Triangle  Park 
_ North  Carolina  27709 


IB.  KEY  WOROS  (Contlnua  on  raaataa  alda  If  nacaaaary  and  Idantlfy  by  block  number) 


National  Science  Foundation 
Washington,  DC  20550 


20.  ABSTRACT  ( Contlnua  on  rararaa  alda  If  nacaaaary  and  Idantlfy  by  block  number) 

Certain  elastic  solids  when  subjected  to  sufficiently  high  loads  abruptly 
change  their  mechanical  properties  and  yet  continue  to  respond  elastically  to 
further  loading.  In  one  dimension  such  mechanically  induced  elastic  phase 
transitions  may  be  due  to  a  non-monotonic  stress-strain  relationship.  This 
appears  to  be  particularly  true  for  certain  mineral  crystals,  such  as  calcite. 

This  work  considers  a  one-dimensional  dynamical  problem  for  a  special  ma¬ 
terial.  The  problem  reduces  to  determining  the  location  of  the  uitemal  moving 
boundary  separating  distinct  elastic  phases.  This  phase  boundary’  is  similar  to 


do,;: 


EDITION  OF  1  NOV  BS  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  flRimi  D«(*  Bntarad) 


a  gas  dynamical  shock  wave.  For  the  problem  considered  here,  this 
phase  boundary  is  shown  to  be  governed  by  a  functional  equation  of 
the  form 

4>(<j>(t) )  +  F(4> ( t) )  +  t  =  o 

for  the  unknown  <f>(t),  whose  F(  )  is  a  known  function  involving 
the  boundary  conditions.  The  unusual  equation  is  derived  by  con¬ 
sidering  the  effect  of  acoustic  waves  repeatedly  reflecting  between 
the  phase  boundary  and  the  external  boundary.  The  equation  is  shown 
to  possess  a  unique  solution  and  is  treated  asymptotically  to  determine 
the  large-time  behavior  of  the  phase  boundary. 


FILMED 

10-85 


DTIC 


